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SUMMARY 


Representative experimental results are presented to show the cur- 
rent status of the panel flutter problem. Results are presented for 
unstiffened rectangular panels and for rectangular panels stiffened by 
1 corrugated backing. Flutter boundaries are established for all types 

of panels when considered on the basis of equivalent isotropic plates. 
The effects of Mach number, differential pressure, and aerodynamic 
heating on panel flutter are discussed. A flutter analysis of ortho- 
tropic panels is presented in the appendix. 


INTRODUCTION 


The problem of panel flutter has become critical with the inception 
of hypersonic and reentry vehicles, particularly when the exposed skin 
panels are not designed to carry appreciable structural loads. For this 
type of exposed surfaces, panel flutter becomes the primary weight- 
establishing design criterion. Some typical skin panels for hypersonic 
and reentry vehicles are shown in figure 1. 

The panel at the upper left of figure 1 is a simple unstiffened 
flat plate. For a multibay region of the exposed surface which is made 
up of sheet and stringers, as shown at the upper right of figure 1, the 
portion of the sheet between stringers is considered to act as an inde- 
pendent panel. For the purpose of this paper these two panel types are 
referred to as unstiffened panels. The two panel types shown at the 
bottom of figure 1 are stiffened by corrugations. The corrugation pitch 
is small and the panel and corrugated backing act as a unit. Both of 
the corrugation-stiffened panels have a much higher stiffness parallel 
to the corrugation than in the other direction. 

The flutter behavior of these types of panels has been the subject 
of many theoretical investigations. Comparison of the predictions of 
theory with experimental flutter data shows marked disagreement for all 
but the simplest configurations. Therefore, experimentally determined 
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flutter boundaries must be relied on for design information. Such 
experimental work is presently in progress and some recent flutter 
results for rectangular panels are presented in this paper. 
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flexural rigidity of equivalent isotropic plate 
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flexural rigidity of orthotropic plate in x-direction 

flexural rigidity of orthotropic plate in y-direction 

twisting stiffness of orthotropic plate relative to x- and 
y-directions 

Young’s modulus of elasticity 
shear modulus of elasticity 
normalized depth of corrugations 

frequency parameter, 



lateral aerodynamic loading 

dimensions of orthotropic plate (see sketch 1 in appendix) 


Z EFF> W EFF 


dimensions of equivalent isotropic plate 
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Mach number 

midplane stresses in x- and y-directions 
integer 

differential static pressure 
dynamic pressure, pU 2 /2 

Rx = / D x 

R x = R x A 2 

R y = N y i 2 /D x 
R y D x 

R = 

y it 2 D xy 
Re 

T 
x o 

t 

^EFF 
U 
W 

w n 


real part 

stagnation temperature 
thickness of material 

thickness of equivalent isotropic plate 

velocity of air flow 
lateral deflection of panel 
deflection function (eq. (3)) 

Cartesian coordinates (see sketch l) 


M 

N X jNy 

n 

Ap 

q 
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mass of plate per unit area 
dynamic-pressure parameter, 2ql^/pD x 
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A cr critical value of A 

v Poisson's ratio 

p mass density of air 

t time 

o> circular frequency 


FLUTTER OF UNSTIFFENED PANELS 


Effect of Panel Aspect Ratio 

Experimental results for the flutter of unstiffened rectangular 
flat panels are shown in figure 2. In this figure the abscissa is the 
panel length-width ratio on a log scale and the ordinate is a modified- 
thickness-ratio parameter which has been shown by theory to be the 
primary panel-flutter parameter. In this thickness-ratio parameter, 
t/l is the ratio of panel thickness to panel length, M is the Mach 
number , q is the dynamic pressure at flutter, and E is the Young's 
modulus of the material. In figure 2 a vertical bar represents the 
test results for panels with a given length-width ratio but with varying 
sizes, materials, Mach number, and temperature. The number by each bar 
is the number of test points and the length of the bar indicates the 
scatter of the results. An envelope curve has been drawn to enclose 
the upper end of the scatter band and hence should be a conservative 
boundary. The flutter region lies below this curve. Note that many 
of the panels tested were very long and narrow; such panels are of con- 
siderable practical importance and flutter data on them have not been 
available until recently. Because flutter amplitudes are usually limited 
by the nonlinear behavior of the panels, panel flutter is considered to 
be primarily a fatigue problem. The severity of the flutter is there- 
fore of interest. In the tests represented by the results in figure 2, 
panels with lower length— width ratios usually exhibited relatively large 
amplitude vibrations and many panels failed before flutter could be 
stopped. On the other hand, the long, narrow panels exhibited mild 
flutter. In some cases these flutter vibrations were difficult to dis- 
tinguish from the background noise without a spectral analysis. 

The thickness-ratio parameter used in figure 2 does not contain all 
the factors that are known to affect the flutter boundaries - hence, the 
scatter. However, this flutter parameter has proven to be very useful 
in comparing the results from a wide variety of test conditions. 
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Effect of Differential Pressure 
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One of the factors not accounted for in the parameter is the effect 
of differential pressure. Previous investigators (see, for example, 
refs. 1 and 2) have shown that small amounts of differential pressure, 
either positive or negative, reduce the thickness required to prevent 
flutter. The results presented in figure 2 were obtained by holding 
the differential pressure as near to zero as possible. This was done 
for two reasons. First, the design information obtained should be con- 
servative, and second, the influence of other factors can be studied 
without being masked by scatter caused by varying differential pressures. 


Effect of Midplane Stresses 

Another important effect on panel flutter is the influence of 
stresses and buckling caused in the panel by aerodynamic heating and 
loading. This effect has been recently investigated in detail by means 
of a series of tests like the one illustrated in figure 3« This figure 
shows the history of a typical tunnel run at M = 3* a stagnation tem- 
perature of 500° F, a dynamic pressure of 3*100 psf, and zero differen- 
tial pressure. During this run the panel temperature varied with time 
as shown. All other test conditions were held constant. The panel did 
not flutter until the temperature reached 150 °, where flutter started as 
indicated by the open symbol. When the panel temperature reached 300° 
the flutter stopped, as indicated by the solid symbol. Theory predicts 
that as a panel is heated the effects of thermal expansion will develop 
compressive midplane stresses that increase the susceptibility to flutter. 
After buckling, the additional temperature rise increases the depth of 
buckles, thereby stiffening the panel and finally stopping the flutter. 

The results of many such runs are shown in figure 4. The thickness- 
ratio parameter for flutter is shown as a function of the panel tempera- 
ture rise. These results were obtained on aluminum panels with a length- 
width ratio of 10 tested at a Mach number of 3* Only panel thickness, 
stagnation temperature, and dynamic pressure were changed from test to 
test. Again the open symbols indicate flutter starts and the solid 
symbols indicate flutter stops. If these results are compared with the 
results shown in figure 2 it can be seen that the maximum value of the 
thickness-ratio parameter for the start of flutter from figure 4 agrees 
with the envelope curve in figure 2 for an z/w of 10. Furthermore, 
the range of scatter of data below the envelope agrees with the range 
obtained by varying the temperature only. Thus, these results indicate 
that midplane stress is the primary cause of scatter. 
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Effect of Mach Number 

The effect of Mach number on panel flutter is shown in figure 5, 
where the thickness-ratio parameter is plotted against Mach number at 
flutter. The tests were run at zero differential pressure on initially 
flat panels buckled by heat. The vertical bars represent the range of 
the test results for each Mach number and the number by each bar repre- 
sents the number of test points. Any substantial Mach number trend in 
these data is obscured by the scatter. Again a large part of this scatter 
is probably due to differences in midplane stress. 


FLUTTER OF CORRUGATION-STIFFENED PANELS 


Next, consider the flutter of the panels stiffened by corrugations 
shown in figure 1. Tests have been run on both types of corrugation- 
stiffened panels shown. The panels were tested both with the flow direc- 
tion along the corrugations and with the flow across the corrugations. 

The flutter results for the corrugation-stiffened panels with the flow 
direction across the corrugations are presented in figure 6. In this 
figure, the thickness-ratio parameter on the ordinate is based on an 
effective thickness, and the length -width ratio on the abscissa is based 
on an effective width. The length l is taken to be the actual length. 

The effective thickness and width are the dimensions of an equivalent 
isotropic plate that would have the same important stiffness character- 
istics as the actual orthotropic, corrugation-stiffened panel. (See 
the appendix.) The equivalent isotropic plate is found to have a greatly 
reduced structural width. Hence, almost -square corrugation-stiffened 
panels have large effective length-width ratios. It should be pointed 
out that this effective length-width ratio is a structural ratio and 
that the aerodynamic length-width ratio of the original panel is unchanged. 

Using these effective dimensions allows the results for the various 
types of panels to be compared on the same basis. In figure 6 the solid 
bars are the flutter results for the corrugation-stiffened panels and 
the shaded bars are the results for unstiffened panels from figure 2. 

The envelope curve from figure 2 has been extended to include the present 
data. This envelope has been drawn as a continuous curve even though the 
results are for entirely different types of panels. At present there is 
no experimental verification that the flutter results for corrugation- 
stiffened panels with effective length-width ratio less than 10 will 
agree with the results for unstiffened panels. The small variation of 
the thickness-ratio parameter over a wide range of the length-width ratio 
indicates that stiffening the panel across the width has very little 
beneficial effect in preventing flutter. Incidentally, some of the 
corrugation-stiffened panels tested represented an actual design for a 
hypersonic vehicle. These panels fluttered at unsatisfactorily low 
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dynamic pressures. A "fix" was obtained on these panels by welding small 
strips to the back of the corrugations, and flutter was thereby prevented 
for dynamic pressures of twice the original flutter value. 

Results for stiffened panels with the flow direction along the 
corrugations do not appear in this figure because no flutter was obtained 
up to the maximum tunnel conditions, which were about 40 percent greater 
than the conditions required for flutter of the best corrugation-stiffened 
panel shown. The panels with corrugations alined with the flow are very 
stiff along the length, and on an equivalent basis the value of the 
thickness-ratio parameter would be very large - well outside the flutter 
region. 

For all the corrugation-stiffened panels that fluttered, the flutter 
was very sudden and very severe. The panels failed unless flutter was 
stopped in a few seconds. 


PANEL FLUTTER ON X-15 AIRPLANE 


The material presented so far summarizes the present understanding 
of the flutter of flat panels based on theoretical calculations and 
experimental results. It is of interest to examine these results in 
terms of a full-scale vehicle. The X-15 was designed as a manned hyper- 
sonic and reentry vehicle and the exposed surfaces of this vehicle con- 
sist of several of the panel types shown in figure 1. Figure 7 shows 
the X-15 and some of the regions where unstiffened panels and corrugation- 
stiffened panels are used. The shaded areas represent fairing panels 
stiffened by corrugations across the flow direction and the black areas 
represent unstiffened panels that have a length -width ratio of 10. Wind- 
tunnel tests on the full-scale vertical tail and on full-scale panels 
from the side fairings have shown that panel flutter exists in the 
operating range of the X-I5. Severe vibrations that appeared during the 
early flights were in part traced to flutter of these side-fairing panels. 

The flutter data for panels from the X-15 are shown in figure 8 in 
relation to the flutter boundary established in figure 6. The ordinate 
is the familiar thickness-ratio parameter based on the panel length and 
effective thickness and the abscissa is the ratio of the panel length 
to the effective width. The circle point represents a flutter point 
for the unstiffened panel on the vertical tail. This was not a boundary 
point (as indicated by the arrow); it was obtained from wind-tunnel tests 
for conditions beyond the design range of the vehicle. The square symbol 
at l /w of 150 was a flutter-boundary point obtained in flight on 

one of the corrugation-stiffened panels from the side fairing, and the 
square symbols at */ w EFF of ^ re P resent flutter -boundary points from 
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wind-tunnel tests of another fairing panel obtained at two different 
Mach numbers. These results are gratifying. The results for the panel 
with ^/ w eff = ^ are seen to agree with the envelope curve established 
in figure 2 by the results for unstiffened panels. This agreement indi- 
cates that the flutter of corrugation-stiffened and unstiffened panels 
can be satisfactorily compared on the basis of equivalent isotropic 
plates. For the panel at ^/ V EFF = a single strap riveted to the back 

of this panel doubled the dynamic pressure required for flutter; three 
straps were necessary to prevent flutter within the design range of the 
vehicle . 


CONCLUDING REMARKS 


Experimental flutter results have been presented for panel configu- 
rations including long narrow panels and corrugation-stiffened panels 
for which previously no information was available. Aerodynamic heating 
was shown to have considerable effect on the flutter results. Flutter 
boundaries were established for all types of panels when considered on 
the basis of equivalent isotropic plates. The results of wind-tunnel 
tests were shown to agree well with flight data obtained from a hyper- 
sonic vehicle. Also, the fact has been established that from the stand- 
point of flutter, corrugation-stiffened panels should be designed so 
that the flow is along the corrugations in order to obtain the greatest 
flutter margin. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., April 12, i 960 . 
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APPENDIX 

FLUTTER ANALYSIS OF ORTHOTROPIC PLATES 


C"— 

O 


i 


The panel configuration to be analyzed herein consists of a simply 
supported flat plate mounted in a rigid wall with supersonic air flowing 
over the top surface. The plate has different elastic properties in the 
flow direction than in the direction across the flow. The plate has a 
length l in the flow direction and a width w and is subjected to 
constant midplane forces of intensities N x and Ny (positive in 

compression) . 



Sketch 1 


In the analysis, small-deflection thin-plate theory for anisotropic 
plates (see ref. 3) is assumed to apply. The equilibrium equation and 
boundary conditions are: 




D x ^ + 2D xy 


d^W d 2 W d 2 W d 2 W , . 

2.2^ — + N *I 2 + Hy ir r ^' L(x ’ y ’ T) 




dx^ dx^dy^ dy* dx^ dy 

W(0,y,r) = W(l,y,T) = W(x,0,t) = W(x,w,t) = 0 


u> 


— (0,y,T) = ll(z,y, t) = ^(x,0,t) = ^K(x,w,t) = 0 
dx 2 dx 2 dy 2 dy 2 

J 

where W(x J ,y,T) is the lateral deflection of the plate and L(x,y,-r) 
is the lateral load per unit area due to aerodynamic pressure. The 
flexural rigidities of the plate in the x- and y-directions are D x 
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and Dy, respectively, D xy is the twisting stiffness relative to the 
x- and y-directions, and 7 is the mass of the plate per unit area. 

It is assumed that the air forces yielded by linearized static 
aerodynamic theory give an adequate approximation. The flutter analysis 
that follows is based on aerodynamic strip theory and follows the proce- 
dure given in reference k for the flutter of isotropic plates. 

The aerodyhamic loading is given by the simple Ackeret value and 
equation (l) becomes 


d iiL + 2d ffL.+ D v + n x + Kv + r &l + £s * , 0 (2) 

X dx 1 * Xy 3x 2 cv 2 y dy 1 * X ' 2 y ■ ‘ 


bxoy c 


dx 2 


ay 2 St 2 P 


where q is the dynamic pressure pU^/2, p = - 1 , and M is the 

Mach number. 


For harmonic motion of the simply supported plate the lateral 
deflection can be written in the form 


W = Re jw n (x) sin e ic0T j 


(3) 


where cjd is the circular frequency. 

Substituting equation (3) into equation (2) yields 

,4 D, 


W n IV - 2n^n 


2 - 2 Qf 5? “n” + + 


- n2 ' 2R y(j) 2H n *»V - - 0 


(<t) 


and the associated boundary conditions 

w n (o) = w n "(o) = w n (i) = w n "(i) = 0 

where the primes denote differentiation with respect to x^ = 


y and 


R v = -5- 
x d y 


N l ^ 

R = 1 

y d v 


A = 


_ 2qP 




PD> 


2 7l 2 
= - cd 


( 5 ) 
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Equation (4) can be rewritten: 

W n IV + AW n " + AW n ' - BW n = 0 


( 6 ) 


where 


A = R x - 
B = k 2 + 


2n 2 * 2 (if $2 
\w J D x 

n 2„2(i\2 R _ n 4 rt 4 £y 

Vw J y Vw ; d x 






J 


(7) 


Equations (6) and ( 7 ) are now in the same form as equations ( 7 ) 
and (8), respectively, of reference k, and the solutions presented 
therein in terms of the parameters A, B, and A apply to the present 
problem. The results presented in reference k- show that the value of 
A cr is essentially dependent on the parameter A and that the results 

for a plate aspect ratio of 1 could be extended to apply to other aspect 
ratios through the use of this parameter. It was shown also that coupling 
between the various modes in the y-direction is slight and that A cr is 

virtually independent of the value of Ry. 


Examination of the parameters of the present problem in equations (5) 
and ( 7 ) shows that the orthotropic plate properties can be replaced by 
the properties of an equivalent isotropic plate that has the same impor- 
tant stiffness characteristics for flutter. From aerodynamic considera- 
tion the equivalent plate must have the same length l and plate stiff- 
ness D x in the flow direction as the orthotropic plate. Since the 
value of A cr is essentially dependent on A, the effective structural 

aspect ratio for the equivalent isotropic plate can be determined from 
the last term in A by means of the relation 



( 8 ) 


where ( — ) is the effective length-width ratio of the equivalent 

W/EFF 

plate . 
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Equation (7) can be rewritten: 


A 


A = R x - 2n^( — 


where 


EFF 


B = k 2 + n c 


(*)La - n \i) 


EFF n 2 
u xy 


( 9 ) 




R x 

R„ = — 
X * 2 


y D xy 


k = 


it 


k_ 

2 


Since A cr has been shown to be essentially independent of 

the essential properties of the equivalent isotropic plate are defined 
by the three quantities Zgjrp, D x , and w^jrp. The last term in B 

appears as a constant and hence would not alter the value of A 


cr 


for 


which the values of B become complex and the motion becomes unstable, 
The properties of the equivalent isotropic plate are as follows: 

l EFF ~ 1 


V EFF 


D- 


EFF 
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= IX 
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xy 
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( 10 ) 




where D 


EFF 


is the plate stiffness of the equivalent plate, t 


is 


the thickness, l 


EFF 


is the length, w 


EFF 


is the width, v 


EFF 

is Poisson 1 s 


ratio, and E is Young T s modulus of elasticity. 
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For the corrugation-stiffened panels discussed in this paper, the 
stiffness properties were obtained as follows: The plate stiffness D x 

was calculated for a unit element of the panel shown in the following 
sketch: 





Sketch 2 

The element consisted of one corrugation and the cover sheet. The stiff- 
ness Dx was obtained by analyzing the element as a structural bent. 
(See, for example, ref. 5*) The stiffness Dxy was obtained from the 
torsional stiffness of the element shown in sketch 2 by means of the 
relation 
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where A Q is the enclosed area of the torque cell, ds is the length 

of a differential element of the sheet, t is the sheet thickness, and 
h is the normalized depth of the corrugation. The properties of the 
equivalent isotropic plate were obtained from equations (10). 
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PANELS FOR HYPERSONIC VEHICLES 
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Figure 1 
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TYPICAL TRANSIENT- HEATING TEST 


M = 3.0 ; Tq = 500°F ; q = 3,100 PSF 



Figure 5 


EFFECT OF AERODYNAMIC HEATING ON 
FLUTTER OF CLAMPED PANELS 
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EFFECTS OF MACH NUMBER 

PANELS BUCKLED BY HEAT; Ap»0 
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FLUTTER OF CORRUGATION -STIFFENED PANELS 
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Figure 6 




REGIONS OF X-15 RESEARCH AIRCRAFT AFFECTED 
BY PANEL FLUTTER 



Figure 7 


X-15 PANEL- FLUTTER RESULTS 
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Figure 8 
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